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ABSTRACT 

Studies of debris disks around white dwarfs (WDs) have focused on infrared 
wavelengths because debris disks are much colder than the star and are believed to 
contribute to the spectrum only at longer wavelengths. Nevertheless these disks 
are made of dust grains which absorb and scatter near-UV and optical photons 
from the WD, leaving a fingerprint that can be used to further constrain disk 
properties. Our goal is to show that it is possible to detect near-UV and optical 
effects of debris disks in the star + disk integrated spectrum. We make theoretical 
calculations and discuss the necessary observational conditions to detect the near- 
UV and optical effects. We show how these effects can be used to infer the disk 
mass, composition, optical depth, and inclination relative to the line of sight. 
If the IR excess is due to a disk, then near-UV and optical effects should be 
observed in only some systems, not all of them, while for dust shells the effects 
should be observed in all systems. 

Subject headings: circumstellar matter - white dwarfs 



1. Introduction 

White dwarfs (WDs) are degenerate stellar nuclei with a mass roughly that of the Sun 
and radii one h undredth that of the Sun; consequently, their surface gravity is ~10 4 greater 



than the Sun's. iLiebertl (119841 ) identified some odd WDs with metal-rich atmospheres. With 
such a powerful gravity pulling the chemical elements toward the stellar nucleus, it is some- 
what unusual to have a metal-rich atmosphere. In fact, the timescales for an element heavier 
than hydrogen or helium to sink are small: ~10 2 yr in WDs with hydrogen atmospheres 



1 Departamento de Fi'sica, Universidade Federal de Santa Catarina, Florianopolis, CP 476, 88010-970, 
Brazil 



- 2 - 



(DAs) and ~10 5 yr with h elium atmospheres (DBs) (IJural 120081 ; I von Hipp el fc Thompson 
20071 IPaquette et al.l[l986h . 



Interstellar material accretion onto the WD surface was one of the first explanations 
for the metal-rich atmospheres. Knowing the diffusion timescales of metals in the stellar 
atmospheres and the metallicity of a given star, it is possible to calculate the necessary 
accretion rate to keep this metallicity constant in time (jKoester fc Wilkenll2006l ). Typical 
values are 10~ 18 to 10~ 15 M yr _1 . These values are too high to be explained exclusively 
by interstellar accretion. Furthermore, if there were accretion from the interstellar medium 
onto the DBs, there should be a lar ge amount of hydrogen pollution in th eir spectra, but 
this pollution has not been detected (IDupuis et al.lll993l ; iFarihi et al.ll2008l ). 



Zuckerman &: Becklinl (119871 ) observed an IR excess in the spectrum of G29-38. The 
shape of this IR excess is a bump which peaks at ~10 ym. Its w i dth is roughly 20 fim 
and can be fitted with a blackbody of T c r ~ 10 3 K. iGraham et al.l (jl990l ) argued that an 



asteroid closely approaching G29-38 could explain this infrared excess. When the asteroid 
orbit reaches the Roche radius it is disrupted and forms a disk around the star. This disk 
is heated up by the stellar radiation and emits in the infrared. The disk material falls down 
continuously on the WD giving rise to the observed metal-rich atmospheres. 

Using a disk model, Jura & collaborators (2003,2007,2008) were able to fit the IR excess 
of many WDs. Through these fits they determined disk physical parameters. With a different 
approach, iReach et al.l (120051 ) fit the spectrum to a thin dust shell model. 



While searching for interacting binary WDs, iGansicke et al.l (120061 . 120071 . 120081 ) observed 
double peaks in calcium lines in a few DAZs and DBZs. Although these observations strongly 
suggest the disk hyp othesis,it is still po ssible that in some cases the emitting region could 
be a torus or a shell (IReach et al.ll2009l ) rather than a disk. 



Previous works focused on the IR emission properties of debris disks around the WD. In 
this work, we propose a new and complementary observational test looking at the absorption 
and scattering properties instead. We develop a simple theoretical framework to predict what 
will be observable and measurable according to the properties of the system. We also suggest 
an observational program to reach our goal. 

We investigate the possibility of detection of debris disks effects in the near-UV and 
optical. We start the analysis in the limit of an optically thick disk in Section [2] and in 
Section [3] we extend this to the optically thin limit case. In Section HJ we discuss the 
observational predictions of our models. Our conclusions are presented in Section [51 
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2. Opaque disk 

The optically thick limit is the natural first approach to investigate the possible effects 
of a debris disk in the spectrum of a WD. This limit can be achieved not only in massive 
disks but also in certain regions of all disks specially if the disk has some gas like some 
recently discovered gaseous disks (IGansicke et al.ll2006l 120071 . 120081 ). Also, the mathematical 
treatment developed in this section will be used in the next section when dealing with the 
optically thin limit. 

A completely opaque disk will not have any spectral features because it is totally opaque 
and absorbs any photon whatever its energy. The only effect of an obscuring disk will be a 
decrease in the received flux from the star. The most the disk can obscure the star is half 
of the projected stellar surface, ttR^/2. The increase in the apparent magnitude of the star 
will be 0.75 mag. 

This value is much higher than current photometric accuracy, and if present, would 
have been detected for those WDs which have an observed parallax, a good spectrum and 
an IR excess. The "good spectrum' permits a determination of T e ff and logg and hence a 
luminosity. Clearly the luminosity inferred from the parallax should agree with the lumi- 
nosity inferred from the spectral fit. At the very least, this procedure will allow us to affirm 
that there is not a big opaque disk in the known WDs with IR excess or, at least, this disk 
is not in a favorable inclination. 

For the general case of a disk with arbitrary inclination and any combination of inner 
and outer radii, the flux received at the Earth from the system (target) is, 



jptarget 



I cos 8dQ 



~D 



tt/2 p<t> 



sin 6 cos 6d(()d6 



where we assumed that the intensity (I) is uniform over the stellar surface. The stellar 
radius is R* and D is the distance from the Earth to the system. There are three possible 
projections for the disk, as seen in Figure [H This gives different values of 4> m in and <p ma x- 



= 27T-2(0 e -&), 
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Fig. 1. — Debris disk and WD. For a given size and inclination of the disk, different amounts 
of it obscure the star. The x— and y— axes are in the sky-plane and the observer is over 
the z— axis which makes an angle i with the normal of the disk. The angles fa and e show 
where the disk starts and stops obscuring the star. 
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, ff™» d<f) = 2 f 7 ' rf0 = 2tt. 

Using the dimensionless radii ?"{i/ e } = R{i/ e }/ Rstar we define the function 5 = rj, r e ): 

7r — (0 e — 0j) , r e cos 2 < sin 9 
g={ ir/2 + 4>i , rjcosi < sin# < r e cosi (2) 

7r , ^ cos i > sin 9 



and write the flux: 



F tar 9 et = j ^ sin ( 2 9) 9 (9, T h T e )d9 (3) 



where 



arctan 



cos i arccos 



sin 2 0/rf e — cos 2 i 



sini 



(4) 



The total flux received from an unobscured system is irI(R*/D) 2 . We call the hypo- 
thetical unobscured star "template" and the obscured star "target" . Defining p as the ratio 
of the obscured to the total projected area: 

-^obscured -^-tarqet / r \ 

P = — A = ~ A — (5) 

Sitotal ^template 

we write the increase in magnitude as: 

Am= -2.51og(l -p) (6) 
In the completely opaque hypothesis we may obviously write that, 



ptarget ptarget / pj 



P = 



ptemplate yj-J y p 



2 



V Jo 



(7) 



sm(29)g(9,r i ,r e )d9 
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The solution of Equation [7] with Equations [2] and [4] gives the flux received from the 
system, as can been seen in Figure [2j The probability of finding a system more inclined 
than a given angle is the ratio of the solid angle occupied by these systems to the total solid 
angle: P(i > i ) = cosi . 

Figure [2] shows that it is possible to detect a completely opaque debris disk. For an 
inclination angle causing any blocking, the bigger the disk is, the easier it is to detect the 
effect. The inner and outer radii are based on physical constrains. If the inne r radius is 
small, the dust particles sublimate because the temperature exceeds ~ 1200K (IJura et al. 



2007al ). On the other hand, if the outer edge of the disk is big (> 100R w d), the dust grains 
will be cold and any emission will be undetectable in practice. The exact values depend on 
the dust type and grain size. 

The presence of an obscuring disk might be inferred by comparing the expected increase 
in stellar magnitude with the luminosity derived from T e g and log g and parallax. If the 
observed and the expected magnitudes are correct and not equal, the flux deficiency can 
probably be explained by obscuration from a debris disk. 

We us e the measured pa rallax of GD 362 to illustrate the previous analysis with one 
real case. iKilic et al.l ( 20081 ) obtained d = 50.6l|j pc for GD 362. Using simple error 
propagation, we have a rough estimate for t he highest a c cepta ble difference between expected 
and measured magnitude: a m ~ 0.15 mag. IKilic et al.l (120081 ) did not find any discrepancies 
between parallax and flux, implying no obscuration of the star. 

From a geometrical perspective this is expected since from Figure [2] Am rs 0.15 mag 
implies an inclination hig her than ~ 8 0° and less than ~20% of the systems will be more 



inclined than this. Indeed, lJura et al.l (j2007bl ) showed that GD 362 must be seen nearly face 
on to be able to reproduce its IR excess flux with physically reasonable inner and outer radii. 
Assuming an almost edge on system would require a big disk and an unusual mechanism to 
heat it to reproduce the measured IR excess flux. Therefore, our work is in agreement with 
the previous results and this analysis illustrates what kind of study must be done with other 
systems which may be found to be nearly edge on. 



3. Optically thin dust disk 

After having derived the basic concepts of the problem with the optically thick limit we 
generalize the equations to the optically thin limit. The ratio of the flux from obscured 
(target) to the equivalent star with no obscuration (template) is composed of three main 
components: 
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Fig. 2. — Increase in magnitude vs. inclination for opaque disks for different inner (r t -) and 
outer (r e ) disk radius combinations. Systems seen face on have % = 0°. The fast decrease 
of Am for i — > 90° is an artifact of the mathematical model that assumes an infinitely thin 
disk. For real disks, there would be a plateau lower than curve peak, but the disk should be 
really thin indeed, so the plateau is very close to Am = and it actually is not a plateau, 
but rather a single point. The top axis label shows the percentage of systems more inclined 
than i. The right axis label shows p, the ratio of the obscured to the total projected area 
(Equation [5]). 
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rptarget 

£ = -La 

^ u ptemplate 

^ u I unobscured ^> u I obscured ^> u I scattered 

The unobscured ratio component is simply 1 — p. The obscured ratio is given by 
p e -T~5 x /cosi fj^g extinction optical depth r^ 2- '* of the disk regions obscuring the star ac- 
counts for the absorbed and scattered light along the line of sight to the star. The scattered 
component comes from the disk regions which do not obscure the star but scatter photons to 
the line of sight. There is no emission component because the dust temperature is lower than 
the dust sublimation temperature (~1200 K) and thus the dust emission only contributes in 
the infrared. 

Using the dimensionless extinction efficiency {Q ext ) instead of extinction cross section 
([C^**] = cm 2 ) we write the differential extinction optical depth in the disk as: 

dr u = nC e u xt dz = nQl xt Tca 2 dz, (9) 

where n (cm -3 ) is the number of dust grains per unit volume, a (cm) is the grain radius and 
z (cm) is the vertical dimension of the disk 

The disk volume density (p) is related to the density of a typical dust grain (p^) through, 

4 

p = -ira 3 p d n. (10) 
Assuming the disk to be vertically uniform we integrate to write 

J-h/2 ^ a Pd ^apd ^apd 

where £ (g/cm 2 ) is the disk surface density and H is the disk height. 
We define: 

3S , . 

ro = — - 12 

and write Equation [8] as 
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& = (1 - P) + pe^V-i + ^| smttered . (13) 

To simplify the scattering term, we assume isotropic and coherent scattering and also 
that the light is not attenuated before and after being scattered by the disk. The last 
hypothesis is valid in the optically thin case and causes an overestimation of the scattering 
because we ignore the absorbed photons. The scattered intensity is given by 

TT TT 

rsca , jwd /-i a\ 

1 = e v : = na Q n.J r , (14) 

COS I COS I 

where e v is the emissivity, Q s u ca is the scattering efficiency and J™ d is the mean stellar intensity 

T wd ir R 2 
Aixr 1 

Ignoring the disk regions hidden by the star we integrate over the disk surface to get 
the flux: 



1 3S ^ arn , (re\ . Twd i R w d 



F "=-2w^yTi) msi * i: (16) 



2 



using Equation [12] and dividing by the template flux, 



UcaUered = lr Qr\n{^)cOS l , (17) 



which allows us to write Equation [13] as 



£„ = (1 -p) +pe -oQ^/co Si + l TQ QSoa ln f™\ cos% (lg) 

2 \riJ 

Besides the parameters p and tq, we have the absorption efficiencies which are char- 
acteristic of the dust t ype. We used the tables of optical constants of silicate glasses from 



Dorschner et al.l (119951 ) . The authors prepared two different glasses in laboratory: pyroxene, 
Mg^Fei-zSiOs, with x=0A, 0.5, 0.6, 0.7, 0.8, 0.95, 1.0 and olivine, Mg2 X Fe 2 -2xSi0 4 , with 
2=0.4 and 0.5. 

Figures [3] and H] display the results from Equation [18] for different optical depths and 
system geometries. Figure [3] represents inclinations where the disk obscures the star, and 
light is absorbed and FigureHJwhen there is no obscuration and we see only scattering plus the 
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WD light. For observational tests, the region from 3000A to 5000A is the most interesting, 
because it shows a sharp change in the ratio between the target and the template which 
cannot be easily discarded as bad flux calibrations. 



Observational test 



Our modeling gives rise to direct observational tests. We can test the effects of the disk 
in the near-UV and the optical, dividing the spectrum of the target by the template. In the 
optically thin case, the result will be color dependent and can provide physical parameters 
for the disk structure. 

The template star should be as similar to the target star as possible. Ideally, it would 
be the same WD without the obscuring disk. As that is not possible to have, we need a 
similar WD without any peculiarity in the spectrum. Any other WD will differ from the 
target star in T e g and log g and this difference can make the division of the spectra resemble 
the expected disk effects. 



In Figure O we present these effects using theoretical WD spectra from lKoester I (120 08). 
We a ssume a target star of T eS = 12, 000 K and log g = 8.0, similar to G29-38 ([Reach et al. 



20091 ). In the upper panel, we keep T eS fixed and var y log.q by 0.05 dex a nd in the lower we 
keep logg fixed and vary T c s by 150 K. According to lLiebert et al.1 (120051 ). the uncertainties 
in temperature are of the order of 1.2% « 150 K and 0.038 in logg. So, larger temperature 
or log g differences would be readily noticed. One can see from Figure [5] that modification 
of UV flux densities by disks that absorb or scatter light, as in Figures [3] and IU can be 
distinguished from observational uncertainties in template stars. 



The White Dwarf Catalog (IMcCook fc Sionlll999l ) currently lists 12,456 stars. Therefore 
it is not too hard to find a template star with a temperature similar to the target. As an 
example, G29-38 and Ross 548 have exactly the same temperature and logg. One needs to 
be careful about this comparison as these values of T e g and log g were obtained from different 
determinations. When comparing the target with template it will be necessary to use T e g 
and logg obtained from similar data and the same models. 



4.1. Parameter determination 

It is difficult to compare directly the definition of p and r with the expected values for 
real disks. We use disk parameters obtained in earlier works to give observational expecta- 
tions and also to help design future observations. 
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Fig. 3. — Expected ratio between the target and the template star (Equation HE]) when the 
inclination is such that the disk obscures a part of the WD as in 1 and 2 in Figure [U The 
disk is composed of olivine Mg .8Fei. 2 SiO 4 dust grains. To calculate the scattering we used 
ln(r e /rj) < ln(100/5) = 3 as a superior limit in Equation [TBI and adopted % = 85°. 
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Fig. 4. — Expected ratio between the target and the template star (Equation IT5|) when the 
disk does not obscure the WD and we see the scattering plus the stellar light. The disk is 
composed of olivine Mgo.sFei^SiC^ dust grains. All the curves were calculated with tq — 0.1 
because the pure scattering term in Equation [IS] is linear in r . 
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Fig. 5. — Comparison of the effects from debris disk obscuration and the effects of small 
differences in T e g and log g between the target and the template stars. The gray lines show 
the expected effects shown in Figures |3] and HI In the upper panel the black solid lines show 
the ratio of WD spectra with fixed T e g- and varying logg. In the lower panel, we fixed logg 
and varied T eff . 
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For disks, the fraction of the obscured to the total projected area, p (Equation EJ), varies 
between when there is no obscuration and 0.5 for a disk which obscures half of the stellar 
surface. However, if the infrared emission region is not a disk but a shell around the star 
( IReach et al.l 120051 ) p will be always 1. Hence, this work provides an independent method to 
test the disk hypothesis. 

In addition to p, we can also determine To- Estimates for the expected values are more 
uncertain, but we can ge t a rough idea by using some mean values for the dust and disk 
properties. iKriigell (120031 ) gives 2.5 g/cm 3 as a typical value for the inte rstellar dust and we 
assume it as a good order of magnitude value for the dust in the disk. lJura et al.l (l2007bl ) 
constrain the disk mass of GD 362 between 10 18 and 10 24 g. Using typical disk sizes of 
10 R w d and 100 R w d for the inner and outer disk radii, respectively, we get a range of r 
from 10~ 4 to somewha t greater than 1, depending on the disk mass and the type of dust 
( iDorschner et al.l 119951 ). Therefore the parameters used in Figure [3] are realistic. 



The disk inclination angle can be inferred from the presence of an flux excess due to 
scattering into the line of sight. Figure H] show s this for inclinations of 0° and 60°. Flux 
excess in near-UV has already been detected by lGansicke et al.l ( 120061 ) in SDSS 1228+1040 
and could be caused by light scattering. For larger inclinations there is a flux deficiency due 
to absorption and scattering out the line of sight, as shown in Figure [3] computed for i = 85°. 
The dividing line between the first or the second case is ~80°. 



5. Conclusions 

In this work, we introduce a new way of looking at the cause of IR excess in white dwarf 
stars. By looking in the near-UV and optical instead of IR we add a new constraint to test 
the disk hypothesis. 

One important distinction of our method, is the fact that the presence of disks would 
cause flux deficiencies in some systems and flux excess in others. We also point out that 
shells would only introduce flux deficiencies effects, and these effects would be detectable in 
all shells. If we find flux deficiencies in every star we observe this would strongly indicate the 
presence of shells rather than disks. Flux deficiencies in only some objects and flux excess 
in others corroborate the idea of a disk. 



If we are convinced disk models are more adequate, d etailed comparisons between disk 
models and data will provide disk mass (I Jura et al.ll2007bl ). composition, optical depth and 
inclination relative to the line of sight. 
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